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ABSTRACT 

We present the cosmological distance errors achievable using the baryon acoustic oscilla- 
tions as a standard ru ler. We begin from a Fisher matrix formalism that is upgraded from 
Seo fc Eisenstein (2003). We isolate the information from the baryonic peaks by excluding dis- 
tance information from other less robust sources. Meanwhile we accommodate the Lagrangian 
displacement distribution into the Fisher matrix calculation to reflect the gradual loss of in- 
formation in scale and in time due to nonlinear growth, nonlinear bias, and nonlinear redshift 
distortions. We then show that we can contract the multi-dimensional Fisher matrix calculations 
into a 2- dimensional or even 1-dimensional formalism with physically motivated approximations. 
We present the resulting fitting formula for the cosmological distance errors from galaxy redshift 
surveys as a function of survey parameters and nonlinearity, which saves us going through the 
12-dimensional Fisher matrix calculations. Finally, we show excellent agreement between the 
distance error estimates from the revised Fisher matrix and the precision on the distance scale 
recovered from N-body simulations. 



Subject headings: large-scale structure of the universe 
ters — cosmic microwave background 



distance scale — cosmological parame- 
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Introduction 



The famous Hubble expansion drives more distant objects to recede faster from us. Recent observations 
of supernovae arg ue that this expansion is in fact accelerating, implying an existence of dark energy with 
negative pressure dltiess et al. 19981 : Perlmutter et al. 1999 : Riess et al. 2001 ; Knop et al. 20031 : Tonrv et al. 
2003t iRiess et al.l 120041 ). This dark energy, which contributes two third of energy density in the present 



Universe, is mysterious in its physical origin. Precise measurements of its time evolution will be crucial to 
uncover the identity of this energy component. One of most promising probe s to measure the dark energ y 
is a standard ruler called baryon acoustic oscillations in large-scale clustering ( Eisenstein et al. 19981 199 



Baryon acoustic oscillations (hereafter BAO) arise from sound waves that propagated in the hot plasma 
of tightly coupled photons and baryons in the early Universe. As the Universe expands and cools, photons 
finally decouple from baryons 400,000 years after the Big Bang, with the sound waves rev ealed as the 



acoustic oscillations in the anisotropies of the Cosmic Microwave Background (hereafter CMB) ( Miller et al 



1999; 
2002: 



de Bernardis et all 2000l: Hananv et al. 2000l: 



Benoit et al. 200 



Bennett et al.l l2003t iPearson et al 



Lee et al 



20011 [Halverson et all 
2003; Hinshaw et al 



2002 



200 



Netterfield et al 



The equivalent 

but attenuated feature exists in the clustering of matter, as baryons fall into dark matter potential well 



after the recombinatio n ( Peebles fc Y u 1970; Sunvaev fe Zeldov ich 1970: Bond fc Efstathi ou 1984: Holt zman 
I — II " II ^' * 'ii i B a 

39c IHu fc S ugivama 1996; Eisenstein & Hu 1998). These were recently detected in galaxy re dshift surveys 

isenstein et alj|200a ICole et alJbood lHutsJl200a IPercival et alJbood : llegmark et alJbooaV Good CMB 
anisotropy data provides the absolute physical scale for these baryonic peaks; comparing this to the observed 
location of the baryonic peaks in the two-point correlation function or power spectrum provides measurements 
of cosmological distance scale. Clustering in the transverse direction probes the angular diameter distance, 
and clustering in the line-of-sight direction probes the Hubble parameter. The cosmological distance scale 
as a function of redshift is then the record of the expansion history of the Universe, which in turn measures 
the evolution of dark energy density. 

Many studies have been devoted to estimations o f precision on cosmological distance scale achiev- 



2003: Hu fc Haiman 


2003: Seo & Eisenstein 2003: Cooravl 12004c iMatsubaral 120041: lAmendola et alT 


2005 


Blake & Bridle! 2005 


Glazebrook & Blake! 2005; 


Dolnev et al. 


2006; 


Zhan & Knoxl2006; 


Blake et al. 


2006) 



A common scheme is to assume the Gaussian statistical errors on the power spectrum, which are constructed 
straightforward fro m the finite volu me of the survey and the incomplete sampling of the underlying density 



field with galaxies (jTegmarkl 119971 ). and then to propagate these errors to constrain the errors on cosmo- 



logical parameters including cosmological distance scale whether by using Monte -Carlo simulations or using 



analytic method such as a Fisher matrix formalism. In ISeo fc Eisenstein! (|2003l) . we used a Fisher matrix 
formalism to show that the BAO from future galaxy redshift surveys can constrain the cosmological distance 
scale to percent level precision, thereby providing robust measurements of present-day dark energy density 
and its time-dependence that are comparable to future Type la supernova surveys. 

The most important element in these methods, other than the survey specifications, is how to include 
the effects of nonlinear structure formation, which depends on redshift as well as bias. The primary non- 



growth of density fields into present-day structures (IV 


eiksin. White. & Peacock 1999: Soringel et al. 


2005 


Aneulo et al.ll2005: 


Seo & Eisensteinll2005 


Whitell2005: 


Crocce & Scoccimarro 


2006: Jeong & Komatsu 


2006 


Huff et al. 20061). This nonlinear mode-coupling also alters the overall shape 


of the power spectrum by in- 



creasing the small-scale power relative to the linear growth of power spectrum (jJain fc Bertschingerlll994 . 
and references t herein). At the same time, this mode-coupling in creases statistical errors above our Gaussian 
approximations ( Meiksin fc Whitell 19991 : Scoccimarro et al . 1999). The nonlinear effects first appear on small 



scales and then proceed to larger scales, tracing the hierarchical formatio n of structures. Nonlinear redshift 



distortions and nonlinear galaxy bias can further degrade the BAO (e.g.. ISpringel et al.ll2005l ; lAngulo et al 



20051; ISeo fc Eisensteinlbooa* Iwhite!l2005l ; iHuff et al-lfeood ) and modify the statistical variance 



In ISeo fc Eisenstein! (|2003f ). we used a sharp wavenumber cut to exclude the nonlinear wavenumbers 
from the Fisher matrix formalism. That is, we treated all wavenumbers up to a certain threshold as linear 
and excluded all larger wavenumbers from our calculations, while the threshold varied depending on redshift. 
Our choice of linear wavenumbers was based on the magnitude of underlying mass density fluctuations rather 
than direct estimations of degradation of baryonic signature with redshift, as we did not have any means 
to estimate the latter without N-body simulations. In addition, the effects of bias or redshift distortions 
on the BAO were too ambiguous to parameterize and include in the Fisher matri x. However, the nonlinea r 
effect does not turn on sharply with wavenumber but instead appears gradually. In I Seo fc Eisenstein! (|2005l ). 
we used N-body simulations to study the effect of nonlinear growth, bias, and redshift distortions on the 
appea rance of the BAO. While we have shown in that paper that our choice of linear scales in Seo fc Eisenstein! 
(|2003l ) reasonably represents the amount of standard ruler information at different redshift, we still lacked 
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a quantitative description of the gradual erasure on the BAO as a function of redshift, bias, and redshift 
distortions. 



Recently, a quantitative model has been provided by lEisenstein et al.l (|2006). In a two-point correlation 
function, the BAO sig nature is an excess of pairs separated by its characteristic scale of ~ 100ft. -1 Mpc. 



Eisenstein et al.l (|2006f ) models the nonlinear process of erasing the BAO signature as di fferential motions 



or Lagrangian displacements of these pairs of tracers. Using the N-body simulations from lSeo fe Eisenstein 
([20051 ). we estimated the amount of differential motions caused by nonlinear growth, halo bias, and redshift 
distortions at different redshifts. We then showed that the nonlinear erasure of the BAO signature is 
successfully modeled by approximating these displacement field as Gaussian and convolving this Lagrangian 
displacement distribution with the BAO peak in the two-point correlation function, or multiplying th e linea r 
power spectrum by the Fourier transform of the Gaussian function. As noted in lEisenstein et al.l ([2006) , 
estimating the Lagrangian displacement fields with a reasonable precision does not require a simulation as 
large as estimating the degradation in the BAO signature does: Lagrangian displacements are calculated 
from all pairs separated by ~ 100ft -1 Mpc, while the BAO depend on the small number of excess pairs at 
that separation. 

In this paper, we use the Lagrangian displaceme nt distribution cal c ulate d from Eisenstein et al. (2006) 
to modify the Fisher matrix formalism presented in Seo fc Eisenstein ( 20031) so as to correctly reflect the 
gradual loss of standard ruler information from the BAO in scale and in time due to not only nonlinear 
growth but also bias and re dshift distortions. We i solate the cosmological distance i nformation from the 
BAO more strictly than in Seo fc Eisensteinl ( 20031 ) . i.e., without information from Alcock fc Paczvnski 
(|1979l ) test or from the broad-band shape of power spectrum, and derive distance error estimates using the 
Fisher matrix ca l culatio ns. We show that the full, 12-dimensional (for one redshift bin) Fisher matrix in 
Seo fc Eisensteinl |2003h can be contracted to a 2-dimensional or even 1-dimcnsional covariance matrix of 
Da and H with physically motivated approximations. We present the resulting fitting formula, for errors 
from the full dimensional Fisher matrix, based on the reduced 2-dimcnsional covariance matrix, which is a 
function of survey volume, galaxy number density, and parameterized nonlinearity. Note that the effect of 
bias and redshift distortions on the BAO is now straightforward to parameterize as an increase in Lagrangian 
displacement field and easy to mingle into the Fisher matrix calculation. We show that the application of 
the formula extends to the photometric redshift surveys. This fitting formula will save considerable amount 
of computational efforts in forecasting distance errors for future galaxy surv eys. We compare t hese results 
with distance constraints f rom a ~x 2 analysis of the full N-body results from I Seo fc Eisensteinl (|2005r) . Our 
fitting formula differs from lBlake et al.l (|2006T ) in many details, such as the different treatments of nonlinear 
degradation of standard ruler test and ours having a structure based on physically motivated model. 

In § [21 we present the full Fisher matrix calculation that is upgraded by using the Lagrangian dis- 
placement field. In § [3J we show a 1-dimensional model as an approximation of the full Fisher matrix in 
§ HI In § SI the 1-dimensional model is extended to a 2-dimensional model. We present the resulting fitting 
formula. We also discuss the extension of this formula for photometric redshift surveys. In § [5j we compare 
distance errors from the fitting formula with the estimates from the full Fisher matrix. In § [6[ we compare 
the distance errors from the full Fisher matrix calculations (or from the fitting formula) with the distance 
errors from a x 2 analysis of the N-body data. 
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Full Fisher matrix calculation 



In this section, we construct the Fisher matrix that assesses the amount of distance information available 
exclusively fr om using the BAO as a standard ruler. This enables us to use the Lagrangian displacement 
distribution (jEisenstein et al. 2006 ) to account for the nonlinear effects on baryonic peaks and therefore 
nonlinear effects on the distance precision. 

Assuming the likelihood fu nction of the band powers of galaxy p ower spectrum to be Gaussian, the 
Fisher matrix is approximately ( Tegmarklll9971 ; ISeo fe Eisensteinll2003 



F — 



V eff (fc); 



dk 



d\nP(k) din P(k) 
~dp~ t KoffW 2(2^)3 

fc max Q\ nP (k^) d\nP{k,ii) 



(1) 



dpt 



dp 3 



V e s(k,fi) 



2irk 2 dkdfj, 
2(2tt) 3 



where P(k) is the observed power spectrum at fc, fi is the cosine of the angle of fc with respect to the line- 
of-sight direction, pi are the cosmological parameters to be constrained, and V e s is the effective volume of 
the survey, given as 



V r eff (fc,/Lt) 



n{r)P{k, ji) 
n(r)P(k,fi) + 1 



dr 



nP(k, /it) 
iP{k,n) + l 



V m 



2\2 



nP(fc)(l + /3/i 2 ) 
nP(k)(l +/3/Li : 



2^2 



f 



v s 



survey • 



(2) 



Here n(f) is the comoving number density of galaxies at the location r, (3 is the linear redshift distortion 
parameter, and Purvey is the survey volume. This formulation assumes that the density field is Gaussian 
and that boundary effects are not important. The second equality in equation f^J) holds only if the co moving 
numb er density n is constant in position, and the third equality assumes linear redshift distortions (jKaisei 
1987h . 



From lEisenstein et al.l ((2006J), the Lagrangian displacement fields due to nonlinear growth and nonlinear 
redshift distortions at separation of lOO/i" 1 Mpc can be approximated as an elliptical Gaussian function. 
Since most of the effect comes from the bulk flow, galaxy bias does not have much effect on Lagrangian 
displacement. The rms radial displacements across (Ej_) and along the line of sight (Ey) at this separation 
follow E^ = EoG and Eii = EoG(l + /), where Eo = 12.4ft, -1 Mpc for a cosmology with <7g = 0.9 at the 
present-day, / = d(ln G) / d(ln a) ~ f2^ 6 , and G is the growth function while G is normalized to be G = 0.758 
at z = such that G(z) = (1 + z) _1 at high z. E should be scaled proportionally to erg. Then the surviving 
baryonic features in the nonlinear power spectrum can be expressed as: 



P b>a l(k,n) 



-Pb,iin(fc,M)exp 



Pfe,lm(fc,/Li) exp 



hi E^ 



fc 2 y 2 ~ 



-k 2 



(l- M 2 )Ei 



,2 v2 " 



(3) 
(4) 



where -Pbjin is th e portion o f the linear power spectrum P\[ n with the acoustic signature, including the effects 
of Silk damping |Silklll968h . 



From equations JT]) and ((2|), the distance precision from a galaxy redshift survey depends on Kj Urve y, 
nP, and the redshift of the survey. The redshift dependence enters because of the loss of information due to 
nonlinear effect, such as the erasu r e of b aryonic features as well as the nonlinear growth of power at large 
wavenumber. In ISeo fc Eisensteinl (|2003l ), the scales degraded by nonlinear effects were removed from the 
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Fisher matrix calculation by cutting off the wavenumber integral at a /c max , where & max depends on redshift 
by requiring er(r) ~ 0.5 at a corresponding r = 7r/2fc max . In reality, the nonlinear effect is progressive with 
k: there is acoustic sca le information even beyond fc max while some information is lost within k < fc max 
( Seo fc Eisenstein 2005h . Now that we have a good model to describe the degradation in baryonic features 



(eq. [3]), we can remove the redshift dependent fc max . We set k max — 0.5h Mpc 1 for all redshifts and change 
S|i and £j_ to account for the erasure of baryonic features. 

As we seek to isolate the information from the acoustic scale, we attach the exponential suppression in 
equation ([3]) to the full power spectrum when computing dP(k, /i)/dpi. This means that we are removing 
all information, baryonic and non-baryonic, from smaller scales, where the broad-band shape of the power 
spectrum might give distance information. Moreover, when we compute the derivatives of P, we take the 
exponential factor in equation f3| outside of the derivatives. This is equivalent to marginalizing over a 
large uncertainty in Sj_ and Ey when computing distance errors; we do not want to include any distance 
information from the anisotropy and scale of the nonlinear damping of the acoustic peaks. 

With these prescriptions, the Fisher matrix becomes 

f 1 f°° dlnP lia (k,n) dlnP hn (k,n) [12^2 ,2 2,^2 v 2 {] ^k 2 dkd^ 

Fij = I J — — Ves(k, /i) exp — fc Sj_ — k [i (T,^ — £_J 



dp, dp, enr,n F L 1 " X/ J 2(2tt) 



(5) 



where Pii n (fc,/z) is the linear power spectrum. 

To this point, we have not excluded the distance infor mation on linear sca l es from non-baryonic features. 



There are two sources of such information: first, from the lAlcock fc Paczvnskil (|1979T ) test (hereafter AP test) 
using angular anisotropy of power spectrum, and second, from the power spectrum not following a simple 
power law. 

First, at a given k, both the redshift distortions from large scale infall and the cosmological distortions 
give rise to angular anisotropics in power. When the linear redshift distortion parameter (3 becomes very 
small, both distortions converge to an identical angular signature, a quadratic in /i, causing the effect of 
redshift distortions to be degenerate from the effect of DaH . In this case, the distance information from the 
AP test will be suppressed. However, when (3 is non-negligible, as we assume a specific angular dependence 
of redshift distortions (e.g., linear redshift distortions), the two effects can be in principle distinguished 
as both distortions have different higher multiples, and we can therefore isolate the cosmological distance 
information on DaH (i.e., the AP test) from the effect of (3. 

Second, if the galaxy power spectrum follows a simple power law without any preferred scale, the 
cosmological distortions remain degenerate from the redshift distortions for a negligible (3. However, the 
galaxy power spectrum contains another preferred scale, the horizon scale at the epoch of matter-radiation 
equality. This e ffect produces much broader feature in power spectrum, but in principle can provide an extra 



standard ruler (|Coorav et al.ll200lh 



However, we consider both of these sources of information less robust than the BAO. For the AP test, we 
do not believe we understand the quasilinear and nonlinear behavior of redshift distortions very accurately. 
For the broadband shape of the power spectrum, the resulting distance information will be susceptible to 
the systematic effect such as tilt, nonlinear bias, or nonlinear growth in power. 

Therefore we want to remove distance information from both non-baryonic features. For the distance 
information from the broad-band shape, we remove this by computing the Fisher matrix -Fh 6= o.oo5 with 
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n b = o.oo£l and subtracting it from the Fisher matrix Fn b with the fiducial f^. In addition, in order to 
remove any distance information from a presumed form of redshift distortions, we assume (3 ~ in the 
derivatives in equation (0 in computing both Fh b= o.oo5 an d the fiducial Fisher matrix. That is, the angular 
AP test is removed. Meanwhile, we hold the amplitude of power in V e g and therefore the statistical error 
per each k bin constant by keeping (3 in the redshift distortions R in V e g (eq. [2]) unchanged. Note that our 
inclusion of an arbitrary growth function already implies that any constraints on j3 are not yielding distance 
scale information from the amplitude of the power spectrum; setting (3 ~ in derivatives will not have any 
further effects. The resulting Fisher matrix after subtraction by Fn b =o.oo5 (while <78,m is the same) will only 
contain the distance information from the BAO. We ignore the minor effect that the broad-band shape of 
the fiducial power spectrum is slightly different from that of the fl b = 0.005 case. 

The resulting Fisher matrix is combined with the CMB information (Planck satellite including polariza- 
tion) and then inverted to give a covariance matrix. We calculate the distance errors marginalized over a total 
of 12 parameters including the angular diameter distance (lnD^) and the Hubble parameter (IniJ) at the 
redshift of the galaxy survey and 10 others: £l m h 2 , fl b h 2 , r, tilt (n s ), T/S, the normalization, lnD^(CMB), 
ln/3, an unknown growth rate (G(z)), and an unknown shot noise. 

In detail, we use the 1st year and the 3rd year Wilkinson Microwa ve Anisotropy Probe (her eafter, 
WMAP1 and WMAP3, respectively) results as our fiducial cosmologies dSnergel et al'j I2003L l2006h . Our 



fiducial model for WMAP1 is then Q m = 0.27, h = 0.72, Q A = 0.73, Vl K = 0, fl b h 2 = 0.0238, r = 0.17, 
n s = 0.99, and T/S = 0. For WMAP3, we use Q m = 0.24, h = 0.73, fl A = 0.76, VL K = 0, Q b h 2 = 0.0223, 
t = 0.09, n s = 0.95, and T/S = 0. 

By adopting S n i (i.e., Y,±_ and S11), we no longer need to specify the redshift of survey, and the distance 
error will only depend on V^rvey, nP, and S„i. We assume a fiducial galaxy survey of V smvey — lh~ 3 Gpc 3 . 
While there is a simple scaling relation for varying Purvey, the effect of nP depends on the value of E n i. We 
try various nP and a large range of S n i- We characterize nP as a value at k = Q.2h Mpc -1 , i.e., nPo.2, where 
P is real-space power. In redshift space, the linear theory redshift distortions enhance the power along the 
line-of-sight direction by R = (1 + Pfi 2 ) 2 against shot noise. 



3. 1-D approximation of Fisher matrix : a centroid problem 

Motivated by the single peak in the correlation function, we next consider whether acoustic distance scale 
precision can be modeled simply as the problem of centroiding a peak in the presence of the noise generated 
from shot noise and the CDM power spectrum. We study this in spherical geometry in this section and 
generalize to the anisotropic case in the next. Nonlinear effects will broaden the peak and therefore increase 
the uncertainty in measuring the location of the peak. We approximate the full Fisher matrix based on how 
well we can centroid the location of the peak, that is, the sound horizon s at the drag epoch when observed 
in the reference cosmology: 



Ft 



In f 



v Bl 



1 



(P(k) 



dP b (k) 



d In s 



4irk 2 dk 
2(2tt) 3 ' 



(6) 



1 We cannot use f2j, = to produce a reasonable CMB information. 
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We first find an appropriate form of dP b (k)/d\ns . If the peak in the correlation function is a delta function 
at r = s a , the Fourier transformation of the peak takes the following form in the power spectrum: 



P b (k) oc 



sin ks 

ks Q 



(7) 



A broadened peak will have the same power spectrum multiplied with an additional damping envelope set 
by the Fourier transformation of the broadened peak shape. For example, if the peak in the correlation 
function is broadened with a Gaussian function due to the Silk damping effect (E s ) and Lagrangian dis- 
placement distribution (S n i), then the convolution in configuration space with the Gaussian function is 
Fourier-transformed to a multiplicative exponential factor in Fourier space. 



P fc (fc )cx !^^exp(--fc 2 E 2 /2) 
ks Q 



(8) 



where £ 2 = E^ 2 + S^. 



In reality, the decoupling between photons and baryons is not instantaneous at the recombination, and 
therefore the resulting Silk-damping effect T>(k) integrated over time deviates from a Gaussian and can be 
better approximated b y exp ( — (fc/ fcsiik) 1 " 4 ) with the Silk-damping scale fc s iik (= 1/S S ) ( Hu fc Sugivamal[l996t 



Eisenstein fc Hulll998f ). We therefore have 



P b (k) oc exp (-(fc/fcsiik) 1 ' 4 ) exp (-fc 2 S 2 ,/2) = exp (-(fcS,) 1 - 4 ) exp (-k*T?j2). (9) 



ks Q 



ks Q 



Starting from physical models of the transfer function gives the same results. Let P b = P — P c 



P{k) = Ak n T 2 (k)=Ak r 



Ak n 



fie 



T c {k) 



^T b (k) + ^T c (k) 



2§^^T b (k)T c (k) 



fit 
S If, 



T b {k) 



(10) 
(11) 



As kT c (k) is a slow function of k relative to T b {k), the leading order term for baryonic features can be 
approximated to 

P b (fc)(x2^^T fc (fc). (12) 

From Eisenstein &; Hu ( 19981 ) (and references therein), T b (k) oc T>(k) sin ks /(ks ) where the Silk damping 
effect T>(k) ~ exp (— (fc/Zcgiik) 1 ' 4 ) = exp (— (fcE s ) 14 ). Including the nonlinear damping returns equation ©. 



P b (k) cx SinfcS ° — ' ^ — ' '~ 2 ^ 2 



ks Q 



•exp(-(/!E s ) L4 )expKE^/2). 



As the amplitude of P b {k) grows with redshift by G 2 , we describe the normalization of P b as: 

sin ks a 



P b (k) = V8vr 2 AoPo.5 



kSn 



■expHfcE^expf-fc'E^) 



(13) 



(14) 



where P0.2 is galaxy power at k — 0.2ft, Mpc 1 at the given redshift. We now differentiate P b (k) to calculate 
dP b (k)/d\ns . 



dP b (k) 
d In s 



= V8^A P . 2 



cos ks 



sin ks 



ks Q 



exp(-(feE,) 1,4 )exp(-* 2 sy2). 



(15) 
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Then the Fisher matrix becomes: 

87T 2 K, 



Fin.. 



' survey 



(P(k) 



-1)2 



AqPq.2 cos ks Q - 



sin ks c 
ks Q 



exp(-2(fcS s ) 1 - 4 )exp (-fc 2 £ nl ) 



4:wk 2 dk 
2(2tt) 3 ' 



(16) 



The comoving sound horizon is ~ lOO/i -1 Mpc and most ofbaryonic information comes from k > 0.05ft Mpc" 1 , 
which keeps ks large over the wavenumber of our interest. Therefore, the sinusoidal terms oscillate rapidly 
relative to the wavenumber dependence of the exponential factor. Treating these oscillations as rapid relative 
to all other wavenumber dependence, we approximate the leading term cos 2 ks as 1/2. The large ks values 
also leave sin 2 ks / (ks ) 2 small relative to cos 2 ks and so we drop this term. These approximations then 
yield 



lo (P(k)/P .2 + (nPo^y 1 ) 2 2 
The resulting error on the location of the baryonic peak is 



-exp(-2(fcS s ) 1 - 4 )exp (-fc 2 S 2 , 



4irk 2 dk 
' 2(2tt) 3 ' 



a* J So = sJF\ 



In S a 



V A 2 

v survey -^0 



dk 



k 2 exp (-2(fcE s ) 1 - 4 ) exp(-fc 2 E 2 1 



\ Po.2 



nPo.: 



-1/2 



(17) 

(18) 
(19) 



The fractional error on the location of the peak from the observed galaxy redshift surveys, ai nSo (= 
er So /s ), is equivalent to the fractional error on the distance estimation when the physical location of the 
peak, that is, the true value of the sound horizon s at the drag epoch, is well known from the CMB. More 
specifically, it is equivalent to the distance information exclusively from baryonic peaks as a standard ruler, 
as we have used only <9P&/<91ns in deriving equation (|17jl . In comparison to the full Fisher matrix (eq. [5]), 
knowing that dP(k)/dlnk in equation §5§ is dominated by dPb(k)/dlnk, we hypothesize that equation (fT7|) 
is a good approximation to equation ([5|) in spherical symmetry. 



4. 2-D approximations of the Fisher matrix 

In this section, we upgrade the spherically symmetric model in the previous section to a 2-D model. In 
allowing anisotropic behavior of the correlation function, we want to measure the location of the baryonic 
peak in the correlation function along (sii) and across the line of sight direction (sj_) in the reference 
cosmology, that is, two axes of an obla te ellipsoidal ridge in the correlation function. The Fourier transform 
of the ellipsoid is sinx/x where x = ^/fc^s^ + k 2 s 2 . 

Measuring the fractional errors on sj_ and s\\ is equivalent to measuring the fractional errors on Da/ s 
and sH, respectively, where s is the true physical value of the sound horizon. When the precision on the 
sound horizon s from the CMB data is much better than the precision on sH and Da/s from galaxy redshift 
surveys, again, the errors on s± and sm are virtually equivalent to the errors on Da and H. 

When Y,± < Eii, the baryonic peak along the line of sight direction is further broadened. In Fourier 
space, the modes along the line of sight are damped further than other modes at given k, introducing an 
angular dependence within the integrand of Fij . 
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We choose parameters p% = lnsj 1 and p 2 = lnsy. We choose lnsj 1 instead of lns^ to be consistent 
with the sign of derivatives for Da in equation l[5]). Using lns^ instead would only change the sign of the 
off-diagonal term of the final Fisher matrix we will derive. The derivatives with respect to the anisotropic 
distances are 



= dlnx/dlnsl 1 = p 2 
f 2 {p) = dlnx/d\ns\\ = p 2 



1 



(20) 
(21) 



when we evaluate the derivatives f\ and f 2 at the fiducial cosmology, i.e., the true cosmology where s± 
sii — s. Using these we can write the 2-D Fisher matrix as 



iv 



v; 



survey 



o (P(k)R((i)+n-iy 



/ 8^A P .2i?(/i) 



OPbao (x) 



d In x 



d In x d In : 

dpi dpj 



2-Kk 2 dk dp, 
2(2tt) 3 



(22) 



where R(p) = (1 + ftp 2 ) 2 is the linear redshift distortions. We can make a similar approximation to equation 
(HSJ for dP b (x)/d\nx: 

(23) 
(24) 



dPb{x)/d\nx oc d(sinx/x)/d\nx 
= cosx — sinx/x. 

We subsequently approximate that cos 2 x ~ 1/2, as before. 
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/p(fc) 

V Po.2 



-k 2 {\- p 2 )T, 2 ± - fcV£ : 



(26) 



where E s is the Silk damping scale, T,± is the real-space nonlinear Lagrangian displacement, and Sii is the 
redshift-space nonlinear displacement along z. Equation (|26[) holds even if the redshift distortion effect R 
deviates from the linear redshift distortions: e.g., R(k,p) (see § 14. ip . The quantity (nP .2-R) _1 becomes 



the effective shot noise N eS . Any additive power due to nonlinear growth or bias will decrease N e g below 
nP Q . 2 R. 



We propose to use equation ([26]) as our fitting function to the full Fisher matrix calculation (eq. [5]). 
As usual, the covariance matrix is the inverse of the Fisher matrix. Unless Ej_ or Em is very large, the p 
dependence of the integrand is mild and therefore easy to compute. We find excellent convergence using a 
simple midpoint method with a grid of 20 points in p and 50 points in k with dk = O.Ol/i Mpc -1 (up to 
k = Q.50h Mpc -1 ). One can precompute P(k)/Po,2 at the gridpoints in k, as these do not depend on nfb.2, 
(3, Ej_, or Eii . The integral over k can be done once for all three matrix elements, so that computing the full 
matrix takes negligibly more time than computing one clement. 

By comparing the numerical results from equation (|26p with those from e quation (l5l), which is pre sented 
in § [5j we derive Aq. We calculate E s (= l/fc s iik) from the equation given in Eisenstein fc Hu (1998 '1- For 



2 fc silk = 1.6(n b h 2 )°- 52 (n m h 2 ) - 73 [1 + (lOAQ m h 2 )- - 95 ] h- 1 ^ Mpc" 1 ) 
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WMAP1, A = 0.4529 and E s = l/k sm = 7.76k- 1 Mpc. For WMAP3, A Q = 0.5817 and E s = l/fc silk = 
8.38ft. -1 Mpc. Note that V SU rvey should be in ft -3 Mpc 3 if k is in ft Mpc -1 units. The derived Aq values are 
consistent with the analytic estimates from equation p0[) although not exact. 

Approximating the Silk damping effect as a Gaussian function, that is, exp (— fc 2 E' s 2 ), is less exact in 
the limit of large E„i and a small nPo.2 but is convenient because we can add the effect of the Silk damping 
and the nonlinear damping quadratically (i.e., E 2 = T,f + E 2 j in eq. [5]). We quote A' Q and T,' s for this 
case as well. For WMAP1, we calculate A' = 0.3051 with a choice of T,' s = 8.3ft -1 Mpc, and for WMAP3, 
A' = 0.3794 with a choice of H' s = 8.86ft -1 Mpc. 



4.1. Redshift distortions and photometric redshift surveys 



The contribution of nP 2 R in equation ([26)1 implies that the effect of redshift distortion R offsets the 
behavior of the shot noise n _1 . On large scales, the anisotropic contribution from P(/i) relative to R = 1 
will increase N e g = uPq^R along the line of sight, decreasing the effective shot noise not only for the 
measurement of H but also for Da, albeit by a smaller amount. 

If the distance to galaxies are uncertain in an uncorrelated way, e.g., due to thermal peculiar velocities 
(i.e., the finger-of-God effect) or due to photometric redshifts, it is straightforward to include this effect in 
equation (|26p by fixing P(/i): we simply need to include a Gaussian uncertainty that corresponds to halo 
velocity dispersion due to the finger-of-God effect or redshift uncertainty for photometric surveys. That is, 

R(k, M) = (1 + /V) 2 exp (-fc 2 /i 2 E 2 ). (27) 

Note that we do not increase values of Em by the amount of uncertainty in the distance. The reason why 
we only need to modify R, whether due to the nonlinear redshift distortions or due to photometric redshift 
errors, is because the resulting exponential suppression in power spectrum (P — > Pexp (— /c 2 /i 2 E 2 )) not 
only decreases the signal but also decreases the variance from the CDM power spectrum in the line-of-sight 
direction. The net effect is thus equivalent to a relative increase of shot noise. In § 15.21 we present distance 
error estimates from a photometric redshift survey using equation (|27j) . 

In the case of photometric redshift errors, there is an additional drawback other than the exponential 
suppression of power in the line-of-sight direction: features in the transverse power spectrum are smeared by 
projections of clust ering at different r edshif ts onto a mean redshift. We ignore this effect in our formulation. 



As we presented in lSeo fc Eisensteinl (|2003h . this effect will increase the error on Da by 13% for photometric 



uncertainty (la) of 4% in 1 + z at z = 1. 

Meanwhile, using P(/i) exp (— fc 2 ^ 2 E 2 ) may be a redundant correction in the case of thermal peculiar 
velocities, as we may already be including these peculiar velocity term in the computation of E|| (of course, 
in redshift space) from simulations. In fact, it is more conservative to put this term into E|| than it is to put 
into E^, because while both cases decrease the signal by the same amount, the latter case also decreases the 
noise by the same amount while the former does not. Therefore, we do not take any steps to isolate thermal 
velocities in the Lagrangian displacements and move that contribution to E z . We simply reserve E s for the 
inclusion of observational uncertainties such as spectroscopic or photometric redshift errors. 
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4.2. A return to spherical symmetry 

When £j_ = Sii and R = 1, the nonlinear exponential damping becomes isotropic, and the 2-D Fisher 
matrix (eq. |26| ) reduces to 1-D Fisher matrix (eq. [17]) multiplied by simple angular integrals. 



F n = F lnSo f {n 2 -l) 2 d[i 
Jo 



(28) 



F22 = F lnSo / {ix'Ydix (29) 



i) 



F12 = F lnSo / -l)d/i (30) 



This gives a 2 x 2 matrix of 







_8_ __2_ 

15 15 



Then the covariance matrix is 



Fa = F )aa A _\ I 5 . (31) 

V 15 15 / 

Cn = = Fl -] o (If). (32) 

The off-diagonal term is not zero, implying that the constraints on sj 1 and sii, and hence the constraints on 
Da and , are not independent even in the limit of superb data. The off-diagonal term being positive means 
that Da and H -1 are anticorrelated. For all modes except those at /i = and fj, = 1, there exist positively 
correlated changes in Da and H that leave the measured quantity of the mode unchanged. Summing 
all the modes, with relative weight depending on /i, leaves a net anticorrelation between Da and if -1 . 
Geometrically, we are attempting to measure the axes of an ellipsoidal shell, but a quadrupolc distortion of 
the shell is less well constrained because it leaves the intermediate angles relatively unchanged. The value 
of r = C12/VC11C22 ~ 0.41 from equation (|32|) . In fact, we find that r w 0.4 regardless of the choice of Sj_, 
Sii , nPo,2, and /3. 

We recommend that this covariance be included in assessing the implications of acoustic scale mea- 
surements for dark energy. We note that Da is an integral of H~ l , so the fact that the acoustic scale 
measurements of these two are anticorrelated means that the constraints are slightly stronger than the 
diagonal errors would imply. 

If one is fitting a model in which the transverse and radial distance scales are required to be the same, 
for example at low redshift, then this implies a contraction of the Fisher matrix with the vector (1,-1). 



This yields an error cr lnQ = cr inDA ^ (1 - r 2 )/(l + 2ra in oJvhxH + of n DA /crf nH ) where of n£ , A = C n and 

a inH — C"22- For the case of spherical symmetry, this reduces to (2/3)a\ n D A (= F^ 1 / 2 ) as expected from 
equation 



5. Testing the approximation 

5.1. Distance errors from the full Fisher matrix and the 2-D model 

We next compare the forecasts on Da and H from the full Fisher matrix (eq. [S], hereafter, 'full-D 
errors') to those from our fitting formula (eq. [26], hereafter, '2-D errors') for various values of nP , 2 , 
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and Sii . We hold R = 1, as this is a minor term. We will present the comparisons for the WMAP3 cosmology; 
the performance is similarly good for the WMAP1 cosmology. 

We first study the case where Eii = Ej_. The upper panels of Figure [1] show expected errors on Da 
and H using the baryonic peaks as a function of nonlinear parameter, Ej_ and En respectively, for various 
nPo.2- The blue lines show the full-D errors and the black lines show the 2-D errors. The figure shows that 
the errors from our fitting formula are in excellent agreement with the errors from the full Fisher matrix 
calculations. Deviations are typically only a few %. The deviation becomes only noticeable for very large 
shot noise nPo.2 (< 0.1), and very large Em (> 20h~ 1 Mpc) where our assumptions to derive equation (f^ 
break down. 

We next present the distance errors when Ej_ < E||. As nonlinear redshift distortions elongate the 
Lagrangian displacement distribution of pairs at 100/i _1 Mpc along z by (1 + /), where / = O ' 6 , the 
distance errors for Ex < E|| represent the effect of nonlinear redshift distortions on baryonic features. We 
consider cases where c = E||/E^ = 2 and 3, while c = 2 corresponds to more realistic redshift distortions, 
and c = 3 depicts an extreme case of anisotropy in Lagrangian displacement fields. 

In the lower panels in Figure [TJ we show the full-D errors and the 2-D errors from equation (|26| for 
c = 2. Again, the 2-D model reproduces the full-D errors to a good extent. We find a similarly good 
agreement between the full-D errors and the 2-D errors for c = 3. 

Note that errors from the fitting formula predict r ~ 0.4 regardless of c and Values of r from the 

full-D calculations are fairly close to the expected values. 



5.2. Photometric redshift surveys 

We use our fitting formula to derive distance errors for photometric redshift surveys. We do this by 
modifying R(n) — > R(k, fi) as described in § 14.11 We use WMAP3 cosmology and assume a photometric 
redshift error (la) of E z = 34ft, -1 Mpc, which corresponds to 1% rms in (1+z) at z = 1. We use equation ([27]) 
to properly include the photometric redshift errors into the Fisher matrix. Figure shows good agreement 
between errors on Da from our fitting formula and errors from the full-D Fisher matrix. Although the 
discrepancy is larger than spectroscopic cases (i.e., Figured]), the offset is at most 8%. Fortunately the 
deviation happens to be small at the common values of nPo.2 for photometric redshift surveys, e.g., nPo.2 = 3 
to 10. 



6. Comparing the distance estimates to N-body data 



In this section, we compare our revised full Fisher matrix formalism (eq. [5]) to the distance estimates 
from a x 2 analysis of N-body data. We find that our revised Fisher matrix formalism provides an excellent 
forecast for the distance errors from the N-body data. 



We use the N-body simulations from lSeo fc Eisenstein ( 2005) and perfo rm a % 2 analysis to extract the 
acoustic scale fr om the simulated power spectrum. Seo &: Eisenstein! ((2005) used WMAP1 cosmology and 
the Hydra code ( Couchman. Thomas, fc Pearcelll995 ) to generate 51 sets of cosmological N-body simulations 
with a box size of 512 3 /i -3 Mpc 3 that were evolved from z = 49 to z = 3 (30 sets at z = 3), z = 1, and 
z = 0.3. In the x 2 analysis, the spherically averaged real-space power spectra Pobs(fcref) from the N-body 
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data are fitted to model power spectra including a scale dilation parameter a, P m (k rc { / a) . In detail, we use 
an additive polynomial function to represent the effect of shot noise and nonlinearity: 



-Pobs(fcref) = (bo + hk xei ) x P m {k lei / a) + (a + aifc rcf + a 2 k 2 c{ ). 



(33) 



The additive polynomial function should also suppress the distance information from the broad-band shape 
of the power spectrum. Recall that, in the Fisher matrix calculations, we explicitly subtracted from the 
Fisher matrix the distance information from the broadband shape (§ [5]). 

Seo fc Eisensteinl |2005h used the linear power spectrum for Pm,(k vc f /a); we improve t his by modifying; 



the lin ear power spectrum to include the nonlinear erasure of the acoustic peaks based on lEisenstein et al 
(|2006h . In detail, we use 



P m = (Plinear(fc) - Psmooth(fc)) exp [-k 2 Y?j2] + P 



smooth 



(*), 



(34) 



where P S mooth is the "no wiggl e" form from Eisenstein fc Hul (j_1998). S n i here does not need to be precisely 
the same as E n i measured in lEisenstein et al.l (|2006T ) and used in the Fisher matrix calculations. S„i is 



used here merely to improve the template power spectrum. The resulting errors on a presented below are a 
smooth function of our choice of £ n i for P m . 



In equation (|3"3"|) , the fit parameters are a, a multiplicative bias bo, a scale-dependent bias b\, and 
additive terms ao, ax, and a 2 from nonlinear growth, bias, or shot noise. The dilation parameter a is the 
ratio of the true distance to our estimate from a \ 2 analysis. For simplicity, we set the true cosmology as the 
reference cosmology, and therefore the mean value of a is expected to be unity. The error on a represents the 
combined errors on -Da(z) and H(z), more specifically a a = (2/3)ctd a (§ I4.2|) , as we use the power spectra 
in real space. 

Since we do not know the true covariance matrix of power spectrum, we do not trust differences in \ 2 to 
give an accurate error on a. Instead, we calculate the mean value and error of a using jackknife subsampling 
of simulations while assuming an error on the band power to be given by a Gaussian random field assumption, 
i.e., based on the number of modes contributing to the band power. We derive 51 (30) subsamples by removing 
one simulation each time from our 51 (30 at z = 3) simulations. Each subsample is then fitted assuming a 
Gaussian random error on the band power. The error and mean value of a is c omputed from the variations 
among the jackknife subsamples. As we emphasized in lSeo fc Eisensteinl (|2005l ). assuming a Gaussian error 
ignores mode-coupled errors between wavenumbers and therefore underestimates the statistical noise on 
small scales. Nevertheless, the variations among jackknife sampling should reflect the true non-Gaussian 
error to a reasonable extent, as these subsamples are drawn from actual nonlinear N-body data. That is, x 2 
statistic slightly misweights the data on small scales in each subsample but the variation among subsamples 
should not produce an overly optimistic a a compared to the true error. 

The resultin g a a from the revised \ 2 analysis is, in general, similar to or slightly smaller than the 
quoted values in Seo fc Eisenstein ( 20051 ). At all redshifts, the value of a a is stable with respect to different 
fitting ranges of wavenu mber k < ka t (&fit = 0.3, 0.4, and 0.5ft Mpc -1 ) or the inclusion of b\. Note that in 
Seo fc Eisensteinl (|2005l ). we had to use a narrower fitting range (fcfit = 0.3ft. Mpc -1 ) at z = 0.3. We will 
quote errors for fcgt = 0.4/i Mpc -1 without b\. When fitting to the matter power spectrum, we will use 
nPo.2 - 80. 

Meanwhile, in all cases but one, the mean values of a are within 1.44cr Q of unity, including 17 out of 
24 cases within \a a , indicating that there is no detection of a bias on the mean value of a. The worst 
out of 24 cases, without b\ with kat = 0.5h Mpc -1 at z — 0.3, gives a mean value of a that is 2.24ctq, off 
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from unity. In detail, in lSeo fc Eisensteinl (|2005T ). we reasoned that using b\ will slightly bias a above unity, 
beyond 1% in worse cases, because the fitting process will favor a negative b\ to match the erased portion 
of the BAO, and the resulting phase shift in the baryonic peaks will be compensated by a above 1. Now 
th at we have accounted fo r the erased features in generating P m (k), this bias on a has decreased compared 
to 



Seo fc Eisensteinl (2005). 



At z = 3, we find a a ~ 0.35% for 



survey 



4ft 3 Gpc 3 from the N-body data. This corresponds to 0.7% 



for Purvey = IK 



3 Gpc 3 . The full-D error (c = 1, S nl = 3.07ft -1 Mpc and nP , 2 ~ 80) predicts a a = 0.65%, 
% s maller than the N-body data. Note that these Fisher matrix errors are different from 



which is about 

the values quoted in lSeo fc Eisensteinl (20051 ). mainly because of the different fiducial cosmology. 



At z = 1, we find a a = 0.38% for V; urV e y = 6.845ft -3 Gpc 



survey 

. = 0.99%. The corresponding full-D error for S n i : 
ii = 5.90 and nPo.2 ~ 80, in excellent agreement. 

At z = 0.3, we calculate <r a = 0.60% for V SI 



For Kmrvey = lh 3 



5.90ft -1 Mpc and nP . 2 



Gpc 3 , this rescales to 
80 predicts a a = 1.01% for 



6.845ft, 3 Gpc 3 of simulations. For V sl 



lft 3 Gpc 3 , the error rescales to a a — 1.57%. In comparison, the full-D error for £„i 
nPo.2 ~ 80 predicts a a — 1.50%. 



8.15ft -1 Mpc and 



We also compute the results when the baryoni c signature in these non linear density fields at z = 0.3 is 
reconstructed using the simple scheme presented in lEisenstein et al.l (|2006h . We use a 10ft _1 Mpc Gaussian 
filter to smooth gravity, displace the real particles and a set of smoothly distributed reference particles by the 
linear theory motion predicted from the nonlinear density fields, find new density fields from the difference 
of the density fields of the real particles and the reference particles, and compute the power spectra of 
these new fields. Fitting these spectra, we find a a of 0.34% for 4ft -3 Gpc 3 of simulation volume^ and thus 
0.68% for lft -3 Gpc 3 . This is a considerable improvement over the 1.57% measured without reconstruction. 
This matches th e Fish er matrix prediction if S„i ~ 3.4ft _1 Mpc, while we measured E n i ~ 4.4ft _1 Mpc in 
( 20061 ). Whether this difference is due to sample variance or due to our technical difficulties 



Eisenstcin ct al 



in estimating E n i for the reconstructed density fields remains to be studied. While this result is an example for 
a negligible shot noise without galaxy bias, our % 2 result implies an impressive prospect for the reconstruction 
in the existing and future galaxy surveys. We will investigate the effects of shot noise and galaxy bias on 
reconstruction in future papers. 

For the biased case at z — 0.3, we use the MASS ca se with m = 10 from Seo fc Eisenstein ( 20051 ). which 
has galaxy bias similar to the LRG samples in SDSS (jZehavi et a.1.1 12005T ) . From the \ 2 analysis, we find 



o a = 0.74% for V sl 



6.845ft -3 Gpc 3 . The error is equivalent to er Q = 1.94% for V sl 



lft- 3 Gpc 3 



To evaluate errors using the Fisher matrix, we need to know the value of nPo.2 that corresponds to this 
case. This is tricky because the nonlinear growth and bias effect in the MASS case with m = 10 increases 
the small-scale power above a nominal shot noise that would be strictly from the inverse number density. In 
addition, this anomalous power is not likely constant in scale. Thus, we roughly estimate N c q = rtPo.2 using 
the best fit additive polynomial in equation p3|l evaluated at k = 0.2ft Mpc -1 : we derive N c g ~ 1.5, which 
corresponds to o a ~ 2.2% for K urve y = 1ft -3 Gpc 3 . 

In general, we find excellent agreement between the errors from the % 2 analysis of the N-body data 
and the analytic full-D/l-D errors. The discrepancy is at most 13% for the cases we have shown. Some 



3 We only include 30 out of 51 sets of simulations because the other 21 sets were produced from initial input power spectra 
that omitted the baryonic signature at fc > 0.3h Mpc -1 . Hence, we use the full 51 only when the non-linearities have reduced 
the role of the k > 0.3h Mpc -1 . 
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of the discrepancy is probably from sample variance. With only 51 values of a, we cannot estimate a a 
perfectly. In addition, while the small-scale power added due to various nonlinearities can be accounted for 
by properly modifying N c g below nPo.2, variance due to nonlincarity are not accounted for in our Fisher 
matrix formalism: this will induce slightly larger distance errors from the N-body data than the full-D/l-D 
errors. While our fitting formula in the x 2 analysis produces reasonable estimates of distance errors, it 



is no t necessarily a mo st optimized one (other examples can be found in iHuff et al.l 12006c iKoehler et al 



2006; 



Smith et al.l 12006). More sophisticated fitting schemes may also reduce the discrepancy between the 



X 2 analysis and the Fisher matrix estimates. 



7. Discussion 

We have computed the cosmological distance errors available from the BAO in future galaxy redshift 
surveys using a Fisher matrix formalism that incorporates the Lagrangian displacement field to account for 
the erasure of the BAO due to nonlinear growth, bias, or redshift distortions. The resulting formalism is only 
a function of survey volume, shot noise, and a nonlinear parameter that can be measured quantitatively. We 
have derived physically motivated lower dimensional approximations to the full Fisher matrix and showed 
excellent agreement between distance error estimates from the approximations and the full Fisher matrix. 
We present the resulting fitting formula to calculate a 2-dimensional covariance matrix for D^(z) and H(z). 
The fitting formula straightforwardly applies to photometric redshift surveys with a simple modification to 
the effective shot noise, although its agreement to the estimates from the full Fisher matrix becomes slightly 
degraded. The merit of the fitting formula is its simplicity in terms of input variables and computation. 

Finally, we compared the error estimates from the revised Fisher matrix with the error estimates from 
a x 2 analysis of N-body simulations. For the x 2 analysis, we also used the Lagrangian displacement field 
to account for the nonlinear effect on baryonic peaks in the template power spectrum. This improved 
various aspects of the error estimates. We showed that both error estimates are in excellent agreement. The 
discrepancies (at most 13%) could be due to non-Gaussianity contributions to the variance not included in 
the Fisher matrix calculation, but may simply be sample variance due to the limited number of N-body 
simulations. 

At http : / / cmb . as . arizona. edu/~eisenste/acousticpeak/baojf orecast .html, we provide a C-program 
for the fitting formula that can be used both for spectroscopic and photometric redshift surveys. 

To use the formula, one must construct Ej_, Em, and nP at the redshift of the survey in question. Ej_ 
and S|| are dominated by the bulk flows in the Universe; their simple scalings wer e given in § [21 although 



one should probably impose a lower <7 8 , matter and rescale E for WMAP3 cosmology. Eisenstein et al. (2006) 



found that highly biased tracers could increase the Lagrangian variances by small amounts, although this has 
not been calibrated in detail. The value of nPo.2 is easy to calculate by using P0.2 = 2710a 2 for WMAP3, 



and P . 2 = 2875cr| . for WMAP1. 



Reconstructing the dens ity field reduces the Lagrangi an displacemen t s and restores the baryonic peaks 
from nonlinear degradation ( Eisenstein et al. 20061 ). In Eisenstein et al. ( 2006h . we found that our simple 



scheme was able to decrease E^ and E|| about by half at z = 0.3 (from Ej_ ~ 8/i _1 Mpc and E|| ~ 
14/i _1 Mpc). In the middle panel of Figure [TJ this decrease due to the reconstruction happens where the 
performance is improving rapidly with E, making future galaxy redshift surveys more promising. The x 2 
results in § [6] implies that the improvements may be better than 50%. We have not yet quantified these 
improvements in E as a detailed function of nP and redshift, but we recommend, to be conservative, that one 
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consider a 50% drop in E as an estimation of what reconstruction can do; Figure 1 of Eisenstein et al. ( 20061) 



estimates that correcting for bulk flow on scales larger than k w 0.1 h Mpc 1 will decrease the displacement 
by 50%. 



Although equation (|26[) appears complicated in that one must compute three two-dimensional integrals, 
it is far simpler than the full Fisher matrix. The 12-dimensional problem requires that one compute many 
integrals of oscillatory integrands to high precision to avoid degeneracies. In contrast, the few integrals in 
the approximation are all smooth and can be computed rapidly and robustly. Meanwhile, the approximation 
includes the anisotropic effects of redshift distortions in the shot noise and the nonlinear degradations of the 
acoustic scale. 

The 1-D model presented in § [3] offers an estimate of how the performance scales with the non-linear 
degradation. Let us rework equation (|19p . replacing the Silk damping form with a Gaussian and approxi- 
mating the denominator as a power-law k 2n . The performance then becomes proportional to 



dk fc 2 - 2n exp(-fc 2 I] 2 ) (35) 

where E 2 = X 2 + E 2 ^. Hence, we find that the distance-scale performance scales as E 1 - 5- ™. For the case 
of white noise, we have n = and the result that the precision scales as E 3 / 2 . This is familiar with white 
noise: a peak has a higher signal-to-noise ratio proportional to the inverse square root of its width, while 
the measurement of its centroid scales as the width divided by the signal-to-noise ratio. 

However, the cold dark matter power spectrum is better approximated by k~ 1 near k ~ 0.15h Mpc -1 . 
Therefore, if our acoustic scale measurements are sample- variance limited, we expect to use n = — 1, which 
implies a precision scaling as £ 5 / 2 . Alternatively stated, the survey volume to reach a given precision scales 
as E 5 . In the language of centroiding a peak, what is happening is that the noise is not white and instead 
has large correlations between neighboring separations. Shifts in the centroid of a narrower peak require 
larger changes on smaller scales that the noise model disfavors. 

This rapid scaling with E implies that one gains rapidly with improvements in reconstruction so long as 
En] is not much smaller than the Silk damping scale of 8/i _1 Mpc. The difference between E n i = 8/i _1 Mpc 
and 4/i~ 1 Mpc is a factor of 1.8 in distance and 3.2 in survey area. It should be noted that the redshift-space 
boosts of the displacements along the line of sight cause E|| to be interestingly large even at z ~ 2. 

Finally, we consider the cosmic variance limits of the acoustic oscillation method. Figure [3] presents the 
fractional errors on Da/s and sH from a survey of 37r steradians, i.e., the reasonably accessible extragalactic 
sky. Redshift bins of width 0.1 are used. The bottom lines show the precision available for a survey 
with perfect reconstruction and no shot noise. The top lines show the precision for a survey with the 
unreconstructed non-linear degradations and shot noise of uPq.2 = 3. The middle lines show a survey with 
reconstruction halving the values of Ej_ and E|| , again with nPo.2 = 3. This last case is what we would suggest 
as a reasonable estimate for a densely sampled survey. One sees that the acoustic oscillation distance scale 
can reach precisions of about 0.4% and 0.7% on Da/s and sH, respectively, for each Az = 0.1 bin at z » 1, 
improving slightly toward z = 3. Of course, most dark energy models predict smooth trends in Da and H 
on scales of Az = 0.1, so in practice one would combine several bins in constructing tests. Hence, effective 
precisions better than 0.2% in distance are available with the acoustic oscillation method. 
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E ± (h" 1 Mpc) 2, (h- 1 Mpc) 

Fig. 1. — 2-D errors from equation (|26p (black lines) and errors from the full Fisher matrix calculations (blue 
lines) for WMAP3. Upper panels: c = = 1. Lower panels : c = 2. Left: distance errors on Da- Right: 

distance errors on H . Off-diagonal terms in the middle field of each panel are defined as r = Cyij \fCv\Gri. 
The bottom field of each panel shows the discrepancy between the 2-D errors and the full-D errors as a 
ratio of the two. The shaded region corresponds to 2% of discrepancy. We find that the 2-D model gives 
excellent fits to the errors from the full Fisher matrix calculations. Solid lines : rtPo.2 = 76.1, short-dashed 
: n Po.2 — 7.61, dotted : nP 02 = 3.81, long-dashed : nP 02 = 0.761, dot-short-dashed : nP 02 — 0.381, 
and dot-long-dashed : nP 2 = 0.076. These values are chosen because nP§ 2 ~ 0.76 is appropriate for the 
Luminous Red Galaxy sample from SDSS in real space. 
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Fig. 2. — Errors on Da for photometric redshift surveys. We assume a redshift error of S z = 34/i _1 Mpc. 
We compare the 2-D errors from equation (f2"6"|) using R(k, fi) in equation (|2"T|) (black lines) and the errors 
from the full Fisher matrix calculations (blue lines) for WMAP3. The lower field of the panel shows the 
discrepancy between the 2-D errors and the full-D errors as a ratio of the two. The 2-D errors from our 
fitting formula are in good agreement with the full-D errors: the discrepancy is at most 8% but smaller for 
nPo.2 = 3 — 10. Solid lines : nPo.2 = 76.1, short-dashed : nPo.2 = 7.61, dotted : nPo.2 = 3.81, long-dashed : 
nPo.2 = 0.761, dot-short-dashed : nPo.2 = 0.381, and dot-long-dashed : nPg.2 = 0.076. 
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Fig. 3. — The fractional errors on Da/ 's {left) and sH (right) available as a function of redshift assuming 
redshift bins Az = 0.1 and a 37r sr survey area. The bottom line in each case is the cosmic variance limit, 
assuming perfect linearity and no shot noise. The top line assumes unreconstructed level of non-linearity and 
a shot noise level of nPo.2 = 3. The middle line uses the same shot noise and assumes that reconstruction 



can halve the values of Si and S 



